DISCRETE ORTHOGONAL POLYNOMIALS AND DIFFERENCE 
EQUATIONS OF SEVERAL VARIABLES 

PLAMEN ILIEV AND YUAN XU 



Abstract. The goal of this work is to characterize all second order difference 
operators of several variables that have discrete orthogonal polynomials as 
eigenfunctions. Under some mild assumptions, we give a complete solution of 
the problem. 
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1. Introduction 

The goal of this study is to characterize all second order difference equations of 
several variables that have discrete orthogonal polynomials as eigenfunctions. More 
precisely, we consider difference operators of the form 

d 

l<ij<d i=l 

where and are the forward and backward operator in the direction of the 
i-th coordinate of M'', respectively, I is the identity operator, A^j , Bi and C are 
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functions of a; G M''. The discrete orthogonal polynomials in our consideration are 
polynomials that are orthogonal with respect to an inner product of the form 

where ^ is a lattice set in M'^ and W is some positive weight function on V. There 
is a close correlation between D, W and V. Some restrictions need to be imposed 
on V due to the complexity of the geometry in higher dimensions. Under some 
mild and, we believe, reasonable assumptions on V, we give a complete solution of 
the problem. 

For d — 1, the one dimensional case, the classification problem was studied by 
several authors early in the last century, we refer to [H [SI for references. It was 
found that the classical discrete orthogonal polynomials, namely, Hahn polynomi- 
als, Meixner polynomials, Krawtchouk polynomials, and Charlier polynomials are 
eigenf unctions of second order difference operators on the real line, and these are 
believed to be the only ones that have such property. To our great surprise, how- 
ever, another family of solutions turns up when we analyze the problem carefully. 
The difference equation satisfied by this family of solutions is similar to that of 
Hahn polynomials but the parameters need to be chosen in a different way. In 
other words, the difference equation has two separate families of solutions when the 
parameters are chosen differently. 

The same phenomenon also appears in the case of two variables. The second 
order difference equations that have orthogonal polynomials as eigenfunctions were 
identified in [12] . but not all families of orthogonal polynomials were found. In fact, 
viewing it as an analogue of second order differential operators that have orthogonal 
polynomials as eigenfunctions, only one family of solutions that had been known in 
the literature was identified for each difference equation. The solutions of difference 
equations, however, turn out to be far richer than that of differential equations. For 
example, in the case of quadratic eigenvalues, only Hahn polynomials on the set 
V = {{x,y) : X > 0,y > 0, X + y < N} are identified in [12]; there are in fact several 
other families, including orthogonal polynomials on Ng, on y = {{x,y) : < x < 
Ni, < y < N2}, and several others. This will be shown in the present paper. 

We start with identifying difference operators D that are self-adjoint with respect 
to the inner product (•,•). This leads to compatibility conditions between the 
coefficients of the operator D. The requirement that the difference equations have 
polynomial solutions and some mild restrictions on V reduce the coefficients of 
D to special simple forms, which allows us to use the compatibility conditions 
to determine the coefficients of the difference operator. A careful analysis of the 
result determines all possible solutions. It is well known that orthogonal bases in 
several variables are not unique [2]. For each solution, we give a family of mutually 
orthogonal polynomials explicitly and compute their norm; in other words, we give 
a family of orthonormal basis explicitly. 

It should be mentioned that the analysis in several variables is by no means a 
straightforward extension of analysis in one or two variables. The complexity of the 
problem increases substantially as the dimension grows. For example, the number 
of solutions grows exponentially with the dimension, and the geometry of the ad- 
missible lattice sets becomes more involved as the dimension increases. Although 
our approach resembles the one used in [12], several new ideas are introduced to 
resolve various outstanding problems. The study in this paper is also much more 
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systematic and thorough, as demonstrated by the new results obtained even in the 
case of d = 2 and by the new family of discrete orthogonal polynomials of one 
variable. 

The paper is organized as follows. The self-adjoint operators are treated in 
Section 2, where the compatibility conditions are derived. In Section 3, we study 
difference equations that have polynomial solutions, called admissible equations, 
and show that the eigenvalues of such equations are either a quadratic polynomial 
or a linear polynomial of the index, and the coefficients of such equations need to 
have certain simple forms under some mild restrictions on V. The compatibility 
conditions can then be used to determine the coefficients, thus D, explicitly. The 
case of one variable is discussed in Section 4, in which the new family of discrete 
orthogonal polynomials is treated in details. For several variables, the case of 
quadratic eigenvalues is studied in Section 5 and the case of linear eigenvalues is 
developed in Section 6. 



2. SeLF-ADJOINT difference OPERATORS 

In this section we discuss the basic properties of self-adjoint second order differ- 
ence operators with respect to an inner product defined on a discrete set. We start 
with some notations. 

Let y be a subset of and let : y ^ K be a positive function. Consider 
the inner product defined by the weight W as follows 

{f,g) = Y.f{x)9{x)W{x), (2.1) 
xev 

for functions f,g:W^^'R. We will be mainly interested in the case when f{x) 
and g{x) are polynomials of a; e M'', but the general discussion in this section does 
not depend essentially on the class of functions. Denote by {ei, 62, . . . , e^} the 
standard basis for R'^. We denote by Ei, Aj and V,, respectively, the customary 
shift, forward and backward difference operators acting on a function f{x) as follows 

Eif{x) = f{x + ei) 

A, fix) = fix + e,) - fix) = iE, - 1) fix) 
^if{x) = fix) - fix - ei) = (X - Er')fix), 

where T is the identity operator. 

We will consider second order difference operators of the form 

d 

D= ^»,jA,Vj +^SiA, + C7I, (2.2) 

l<i,j<d i=l 

where Aij , Bi and C are some functions of x. Immediately from the definition one 
can see that D can also be rewritten as 

d d 

D= aijEiEr^ +J2PiEi + Y,7iE-^ +SI, (2.3) 
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where the new coefRcients aij, /3i, 74 and S are related to the old ones via the 
formulas 

auj = -Ai^j for I j <d (2.4) 

d 

f}^^Yl ^'.fc + B^^oyl<^<d (2.5) 
fc=i 

d 

7, =^Afc,i (2.6) 

k=l 

d 

6 = C- J2 (2-7) 

l<i^j<d i=l 

The representation (j2.2p will be more convenient when we deal with polynomials, 
because the operators A.; and Vi decrease the total degree of a polynomial by one. 
However, necessary and sufficient conditions for an operator D to be self-adjoint 
with respect to the inner product ()2.ip are much simpler and natural if we write 
it as in (|2.3p . Formulas (|2.4p - (|2.7p allow us to go easily from one representation to 
another. 

Let us define "directional" boundaries of V as follows: 

d^V ^ {x eV ■.x±ej iV} for j 1, 2, . . . , d 

d.^jV = {x : X + e,- Ej (^V} for 1 < i ^ j < d. 

Example 2.1. Let No denote the set of all nonnegative integers. For the set 
= N[J we have O'V = d^-jV ^ V n {x : Xj ^ 0} and d+V = for j 1, 2, . . . , d 
and i ^ j- 

Other sets that play crucial role later are listed below 

Vf = {x ^n"}, : x^ <h for i ^ 1,2, . . . , d}, (2.8a) 

where U are positive integers; 

1/^ = {x e NfJ : xi + ^2 + • • ■ + < N], (2.8b) 
where A'^ is a positive integer; 

V^^S ^V^{^{x:x<k for i £ 5}, (2.8c) 

where S' is a nonempty subset of {1, 2, . . . , d} and U are integers such that 1 < < 
N. 

For example, for the parallelepiped we have djVf" — Vf" r\ {x : Xj = 0}, 
d+V^'^ ^V^'^nix: X, = I,} and a,,,V,<^ = d-V^'' U d+V^'^. 

For defined by (p^Sb)) we have d^^V^ ^ ^^JV^ = n {x : Xj = 0}, and 
9^^^ = Vi''n{x:xi+X2--- + Xd^N}. 

The next proposition characterizes the self-adjoint operators with respect to the 
inner product (|2.ip in terms of their coefficients. 
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Proposition 2.2. The operator D is self-adjoint with respect to the inner product 
()2.H) . if and only if 

W{x)-f^{x) ^W{x~e,)(3,{x-e^), x,x-e,eV (2.9) 

W{x - ei)aij{x - Ci) = W{x - ej)aj^i{x - ej), x -e^^x- Cj e V, j ^ i, 

(2.10) 

and 

7j(a;) = 0, xedjV, (2.11a) 
/3j(a;) = 0, xed+V, (2.11b) 
a,j(a;) = 0, x d^,jV, i^j. (2.11c) 

Proof. Using {x-\- ei — Cj : x €V\ di^jV} dj^iV and changing tlie summation 

index we see ttiat 

{aijEiE~^u,v) ^ ^ aij{x)u{x + ei - ej)v{x)W{x) 

+ aij{x + ej — ei)u{x)v{x + Cj — ei)W{x + Cj — ei). 

xev\^J,^v 

Writing similar equalities for {PiEiU,v) and ('-fiE^^u,v), using {x ± Ci : x £ V \ 
dfv} = V\d^V, and comparing {Du, v) with (u, Dv) gives the stated conditions. 

□ 

As an immediate corollary from the relations established above we can deduce 
compatibility conditions that need to be satisfied by the coefficients of the operator 
D. 

Corollary 2.3 (Compatibility conditions). // the operator D is self-adjoint with 
respect to the inner product defined by (|2.ip then for 1 < i ^ j < d we have 

ai,j{x - ei)(3j{x - ej)j,{x) = aj,i{x - ej)Pt{x - 61)7^(2;), (2.12a) 
for x,x — ei,x — Cj e V and 

I3i{x - ei)f3j{x - - ej)7j(a;)7^(a; - e^) = 

(3j{x - ej)(3i{x - Ci- ej)j,{x)"fj{x - e,), (2.12b) 
for x,x — ei,x — ej,x — Ci — Cj £V . 
Proof. Writing 

W{x - e-j) _ W{x - Cj) W{x) 

W{x - Ci) ~ W{x) W{x - e,) 
and using (|2.10p on the left side and (|2.9p for the right side we obtain (|2.12ap . 
Similarly if we use (|2.9p for all ratios in the identity 

W{x) W{x-e^) _ W{x) W{x-ej) 
W{x - Ci) W{x - Ci - Cj) W{x — Bj) W{x - Ci- Cj) ' 
we get (|2.12bp . which completes the proof. □ 

Remark 2.4. The compatibility conditions stated above can be easily extended 
to more general difference operators. For example, we can consider the operator 
D that adds two additional terms EiEj and E~^E~^ in the operator D in (|2.3p . 
However, the complexity of the computations in the sections below will increase 
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significantly even with these two terms added. Furthermore, since A = V + VA, 
the operator D includes the 4-th order term A1V1A2V2 (say, d — 2) which is the 
discrete analog of where di stands for the partial derivative with respect to 

Xi . Classifying orthogonal polynomials that are eigenfunctions of D appears to be 
an interesting but much harder problem. 

3. Admissible equations and orthogonal polynomials 

In this section we study difference operators D that have discrete orthogonal 
polynomials as eigenfunctions. For more details on the general theory of discrete 
orthogonal polynomials of several variables we refer the reader to [TT] . 

Throughout the paper we use the standard multi-index notation. A multi-index 
will be denoted by /i = (^1,^2, ■ • ■ G Nq- For each ^ we denote by a;^ the 
monomial 

— \ 2 d 

of total degree = /-ii + fJ.2 + ■ ■ ■ + fJ-d- The degree of a polynomial is defined 
as the highest degree of its monomials. We denote by R[x] = M.[xi,X2, ■ ■ ■ , xj] the 
space of all polynomials in the variables xi,X2, ■ ■ ■ ,Xd and by Hf^ the subspace of 
polynomials of degree at most n in the variables {xi,X2, ■ ■ ■ ,Xd}- The latter has 
dimension dim(n^) = ("+''). 

Let V be an at most countable set of isolated points in R'' and let \V\ denote the 
cardinality of V. If (/, g) — for the inner product defined by p.ip . we say that / 
and g are orthogonal with respect to W on the set V. Orthogonal polynomials on 
V depend on the structure of the ideal 

3{V) — {p £ R[a;i, a;2, . . . ,Xd] ■ p{x) — for every x £ V}. 

In fact, the orthogonal polynomials belong to the space 

R[V]^R[xuX2,...,Xd]p{V). 

There is a lattice set A{V) such that every polynomial in M.[V] can be written as 

AiGA(V) 

The lattice set is not uniquely determined by V. However, in all cases that arise 
here, there is a natural way to identify V and A{V), see Example 13.11 Denote 
Afe(F) = {/i G A(F) : = k} and let r^. be the number of the elements in the set 
Afc(F). 

Example 3.1. It is clear that for 1/ = Nq we have 3(Nq) = (0) and therefore 
R[Nq] = R[a;]. For all other sets defined in Example 12. II one can easily find (induc- 
tively on d and \V\) explicit sets of generators of the ideal J(V^). In particular, the 
generators listed below allow us to identify A{V) and V. 

For VJ"*, defined by (|2.8ap . the ideal 3(V;'') is generated by the set 

Gf = {{-^1)11 + 1, i~X2)l2 + l, . . . , {-Xd)l^ + l}. 

For V^, given by (j2.8bp . the ideal 3(V^) is generated by the set 
Gff = {(-xi)pi (-X2)^, • • • i-Xd)^, : ImI - ^ + !}• 
Finally, for g in (|2.8cp , the corresponding ideal is generated by 

G%^s = G%U{{~x,)u+i:teS}. 
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Next we want to consider difference operators D on the space M.[V]. In order to 
have a proper action of D on M.[V] the ideal 3(V) must be I?-invariant, i.e. 

DiJiV)) c 3{V). 

One can easily show that if D is self-adjoint with respect to the inner product 
()2.ip . then the above condition is satisfied. More precisely, we have the following 
proposition. 

Proposition 3.2. If the operator D, written in the form (|2.3|) . satisfies conditions 
(|2.1ip then the ideal 3{V) is D-invariant. 

Proof. Notice that if a polynomial p{x) vanishes on V, then the polynomial p{x — ej) 
wiU vanish on 1/ \ d~V. Thus, if p{x) e 3{V) and (|2.11ap holds then ^j{x)p{x - 
Cj) G '3{V), i.e. the operator ^jEJ^ preserves '^{V). Similar arguments show that 
(j2.11bp (resp. p.llcp ) implies that (3jEj (resp. ai^jEiEJ^) preserves 3{V), which 
completes the proof. □ 

Definition 3.3. Let 3(F) be _D-invariant. The equation 

Du = Xu (3.1) 

is called admissible on V if for any k ^ Nq there is a number At such that the 
equation Du = X^u has linearly independent polynomial solutions in R[V"] and 
it has no nontrivial solutions in the set of polynomials of degree less than k. 

Remark 3.4. Notice that if ()3.ip is admissible on V and if we have nontrivial 
solutions of the equations Du — XkU and Du — Xiu for distinct integers k,l e Nq 
then Xk ^ Xi. Thus if D is also a self-adjoint operator with respect to the inner 
product (j2.ip . the polynomial solutions for Du = XkU and Du = Xiu will be 
mutually orthogonal. 

From now on, we will assume that Aq = (which is equivalent to C = for 
operators D of the form (|2.2p ). This is not a real restriction on D since one can 
replace D by D — AqI if necessary. 

The next proposition shows that we can pick linearly independent polynomial 
solutions in such a way that the highest terms contain single monomials. 

Proposition 3.5. The equation p.ip is admissible on V if and only if for each 
fc e No there exists a number Xk such that the equation Du = XkU has r^ linearly 
independent polynomial solutions of the form 

P^{x)^x'' modn^_i. (3.2) 

The proof of Proposition 13.51 follows easily form Definition 13.31 see [T^ in the 
case d = 2. 

Remark 3.6. If D is an admissible operator on V we can trivially extend it to an 
admissible operator on the set V — {{x,0) : x ^ V} C K'^+'^. When we classify the 
possible admissible operators, we will exclude such degenerate situations. 

The next proposition essentially characterizes the admissible difference operators 
of the form (|2.2p . In order to eliminate several singular cases we will assume that 
fj. E A(F) for 1^1 = 1,2,3. Instead of putting fi e A(y) for |^| = 3 we can require 
that 3ei G A(y) for some i G {1,2, ... ,d} and for all 1 < i 7^ j < d at least one 
of {2ei + Cj, Ci + 2ej} belongs to A{V), but this would make the statement of the 
theorem long and awkward. 
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Theorem 3.7. Assume that /i G A(V") for = 1,2,3. Let D be a second or- 
der difference operator defined by ()2.2p such that 3{V) is D-invariant. Then the 
following conditions are equivalent: 

(i) The equation Du = Xu is admissible; 

(ii) The coefficients Bi{x) are polynomials of degree at most 1, and Aij(x) are 
polynomials of degree at most 2 satisfying 

Bi — bxi mod IIq for i = 1, 2, . . . , d (3.3a) 

A,^, = axi mod /or i = 1,2, .. .,d (3.3b) 

Ai,j + Aj^i — 2aXiXj mod Ilf for 1 < i ^ j < d (3.3c) 

for some constants a and b. The eigenvalues Xk must be distinct and they 
are given by 

Afe = k{ka -a + b). (3.4) 

Proof. Assume first that the equation Du — Xu is admissible. Applying Proposition 
13.51 we see that there exist polynomials 

Pe . {x) = Xi mod IIo 

satisfying 

Du = Xiu. (3.5) 

Notice that A^-Pe. = and AjVuPe, = for all j,k = 1,2, ... ,d. Thus DPe, = 
Bi(x) and the last equation shows that formula (|3.3a[) holds, where we put 5 = Ai. 
Similarly, we can find polynomials of the form 

P2eAx)^xj modn^ (3.6) 

satisfying 

Du = X2U. (3.7) 

This time we have AjP2ei — 2xi6ij mod IIq and A^VfePje^ — 26ij6i^k- Using now 
(I3.3ap we see that 

DP2e, = 2A,,j + 2bxf mod Uf, 
which, combined with (|3.6p and p.7p . shows that 

Aut = ^ ^^ a:- mod H^, 

thus proving (j3.3b[) if we denote a = (A2 — 2b) /2. 
Next, for i =/= j we take a polynomial of the form 

Pe-+ejix) = XiXj mod uf 

satisfying (|3.7p . A simple computation as above gives that 

DPe^+e^ = Ai,j + Aj_i + 2bxiXj mod Hf, 

which combined with p.7p shows that (|3.3cp holds. It remains to show that Aij 
are polynomials of degree at most 2 and that (|3.4p holds. First we use a polynomial 
of the form 

Pge, = xf mod 

for some i satisfying 

Du = X3U, (3.8) 
to show that p.4p holds for fc = 3. Indeed, we already know that if P is a polynomial 
of degree at most 2, -D(P) will also be a polynomial of degree at most 2 (because 
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we can write it as a linear combination of 1, Pe^ , Pe^+ej , for i, j = 1, 2, . . . , d). Thus 
D{P3ei) = D{x^) mod nf. On the other hand 

= 6xiA^^, + {3x^, + 3a;, + = (6a + 35)a;f mod H^, 

where in the last equality we used (|3.3a|) and (|3.3b[l . Comparing now the coefficients 
of x^ on both sides in p.Sp , we see that A3 = 6a + 36, which is exactly p.4p for 
A: = 3. 

Next we consider the solution to equation (|3.8[) of the form 

P2e,+ej (x) = xjxj mod 

for i 7^ j. Again Z3(P2ei+ej ) = D{xfxj) mod 112, but this time 

D{xjxj) = {2xi + l)Ai^j + {2xi - l)Aj^i + 2xjAis + [2x^ + l)xjBi + xfBj 

= 2xi{Ai^j + Aj^i) + Ajj - Aj,i + 2xjA^^^ + 2xiXjBi + x^Pj mod li^ 

= (6a + ih)xlxj + Aij - Aj.i mod IIj 

= X'ixjxj + Ai.j — i mod IIj, 

upon using (|3.3p and A3 = 6a + 36. Equation (|3.8p shows that Ai,j — Aj^i — 
mod 112, which combined with (|3.3cp proves that Aij and Aj^i are polynomials of 
degree at most 2. Finally, let /i G A(y) such that |^| = k and let 

P^ = xf" mod n^_i 

be a solution to 

Du = XkU. (3.9) 
The admissibility implies that D{u) — D{x^) mod Notice that 

A.x^ = Mi^:^"^' mod n^_2 for all i (3.10a) 
A.Vjx'' = ^ii^ljx''''''-''' mod n^_3 for i 7^ j (3.10b) 
A.ViO;^ = Aii(Aii - l)x''~2^' mod H^.g for aU i (3.10c) 

and therefore 

d 

l<i<j<d 1=1 
d 

+ Y /^^a;'""^* Bi mod U'l_^ 



=1 

d 



(3.11) 



2a ^ ^i^j +a^^i(/ii - 1) + 6^/ii mod nj?^i 

l<j<j"<rf i=l i=l J 

= |^|(a|^| - a + 6)a;^ mod n^.^, 

which combined with p.9p gives (|3.4p 

Conversely, assume that (ii) holds. One can easily show by induction on that 
for fi there exist polynomials of the form 

P^{x) = x^'+ Y Ct^,Mx) (3.12) 
kl<lM 
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satisfying ()3.9|) with eigenvalue A|^| given by ()3.4|) . Indeed, let us assume that this 
is true for \iJ.\ < k — 1 and take fj. such that = k. Using computation (|3.11|) we 
see that 

kl<lM 

for some constants 7^,1/. Then 

\i^\<\fj.\ 

i.e. if we pick Cf_i^i/ — — TM-i^ 

/(A|^l — A|^|) the polynomial defined by (|3.12p will 
satisfy (|3.9p . Thus, in the quotient space M.[V] we have rfc polynomials of total 
degree k satisfying (|3.9p . □ 

In the next sections we will use conditions (|3.3p and the compatibility conditions 
p. lip and (|2.12p to determine the possible self-adjoint operators on appropriate 
sets V. 

Since the difference equation Du = Xu is invariant under translations x — > x + 
we can consider the difference equations modulo translations. 

In the rest of the paper we will impose certain boundary conditions on the dif- 
ference operator D. More precisely, we will assume that for every j the coefficients 
Aij{x) vanish on the hyperplane {x : Xj = 0}. Roughly speaking, this means that 
V has sufficiently many boundary points belonging to {a; : xj — 0}. We show 
below, that if a 7^ and if V has "enough" points on its boundary then, up to a 
translation, this must be true. 

Let us denote by djV — djV r\i^j dijV the "lower" jth boundary of V. We say 
that V has a nontrivial boundary if for every j = 1,2, ... ,c?, djV contains more 
than 2'* points and that these points do not belong to a variety of dimension less 
than d ~ I. For example, a line in rf = 2 is a nontrivial boundary, so is a plane in 
d = 3. 

Proposition 3.8. Let V be a discrete set with a nontrivial boundary and let D be 
an admissible, self-adjoint operator. Ifoy^O, then after an appropriate translation, 
djV becomes a subset of {x : Xj = 0} and the coefficients Aij{x) vanish on {x : 
Xj = 0}. 

Proof. Using the fact that D is self-adjoint. Proposition 12. 2i (|2.4p and (|2.6p we see 
that 

(x) =0, xedjV for l<j<d. 
Applying Bezout's theorem for polynomials in d variables [7, Chapt. IV, Section 
2] for the polynomials Aij{x), . . . ,Adj{x) which satisfy (j3.3bp - p.3cp and vanish 
on djV, one concludes that they must contain a common linear factor. Equation 
(I3.3bp shows that this factor must be of the form Xj — hj. Thus, applying the 
translation x —^ x + h with h = [hi, h2, . . . , hd), Aij takes the form 

Aij{x) = Xj X (linear polynomial in xi,X2, ■ • ■ , Xd), (3.13) 
and completes the proof. □ 

The compatibility conditions also imply that Aij — Aj_i G Ilf for i ^ j, which 
combined with (|3.3cp shows that Ai^j ~ axiXj mod 11!^ (the argument works si- 
multaneously for a ^ and a = 0). The form of the coefficients determined below 
will be the staring point in the next sections. 
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Proposition 3.9. Let Aij and Bt satisfy the conditions in Theorem \ 3. 7]f ii ) and 
assume that Aij vanishes on {x : xj = 0}. Then, the polynomial identities (I2.12ap 
imply that 

Ai,j = Xjiaxi + li_j) (3-14) 

B,, = bxi + Sj, (3.15) 

or equivalently: 

ai,j = -Xj{axi + kj) (3.16) 

d d 



I3i^ axi'^^Xk + '^^li,kXk'Vbxi + Si (3.17) 

k=l k=l 
d d 

7i = axi^Xfe + (3.18) 



for some constants a, h, lij and Si. 

Proof. Equation (I3.15P follows immediately from 113. 3ap . From equations (|3.3bp . 
p.3cp and the fact that Aij vanishes on {x : Xj = 0} we deduce that 

Aij = Xj{qi^jXi + Zij), (3.19) 

where i = a and Qij + qj i — 2a. It remains to show that Qi j = a for all i ^ j. 
From dOIl-dll]) we see that 

aij = —Qi^XiXj + linear terms 

n 

Pi = Xi qi,kXk + linear terms 
fc=i 

n 

7i = Xi Qk.iXk + linear terms. 

k=l 

Comparing the coefficient of x^x^ on both sides of (|2.12a[) we obtain 

which combined with g^.j + qj^i — 2a gives qij — qj,i — a. □ 

In the rest of the paper, we will find the general solution of the compatibility 
conditions (|2.12p . by equating the coefficients of the different powers of x, and we 
will determine the possible sets V, weights W as well as explicit bases of orthogonal 
polynomials. 

Notice that (13. 4p implies that Xk must be at most quadratic in k. Thus, we have 
essentially two possible cases: 

(i) Xk is quadratic in k (i.e. a ^ 0); 

(ii) Xk is linear in k (i.e. a — 0). 

Since we can always divide the equation Du = AfcU by a nonzero constant, we can 
take a = — 1 in (i) and similarly we can normalize so that 6 = — 1 in (ii). These two 
cases are discussed in Sections [5] and [6] respectively. In the next section we treat in 
details the case d—1. 
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4. One dimensional orthogonal polynomials 

The one dimensional case is the simplest case with no compatibility conditions. 
Our goal in this section is to provide a detailed study in this simplest situation and 
to show how the procedure works (i.e. how we determine the weight function W{x) 
and the set V) . The discrete orthogonal polynomials of one variable will also serve 
as building blocks in the higher-dimensional cases. 

Historically, orthogonal polynomials satisfying difference equations of one vari- 
able were studied early in last century, we refer to [TJ [5] for references. It is well 
known that there are four families of such orthogonal polynomials of one variable, 
namely, Hahn polynomials, Meixner polynomials, Krawtchouk polynomials, and 
Charlier polynomials. However, to our great surprise, another family of orthogonal 
polynomials of discrete variable shows up in our study. 

For simplicity of notation, we will drop all indices within this section, i.e. x = 
(xi) G R, E = El, A = Ai etc. The operator D will be of the form 

D = A{x)AV + B{x)A = f3{x)E - {P{x) + A{x))I + A{x)E-\ (4.1) 

where 

A{x) = x{ax + I), B{x) =hx + s, and (3{x) = A{x) + B{x). (4.2) 
As we explained at the end of the previous section, we have two possible cases: 

(i) quadratic eigenvalue A^, we can fix a = —1 

(ii) linear eigenvalue A^, we can fix a = 0, & = —1. 

4.1. Quadratic eigenvalue: a = — 1. In order to make the formulas more sym- 
metric, let us write A{x), (3{x) and B{x) as 

A(x) = -x{x + as) (4.3a) 

P{x) ^ -{x + ai + l){x + a2 + l) (4.3b) 

B{x) = P{x) - A{x) ^ (as -ai-a2- 2)x - (1 + ai)(l + aa), (4.3c) 

i.e. we have replaced I — —03, b — — ai ~ a2 ~ 2 and s = —(1 + ai)(l + 02)- 
Thus (HH) gives 

W{x) ^ (i{x - 1) ^ {x + ai){x + a2) 
W{x-l) A{x) x{x + a:i) ' ^ ' ' 

There are essentially 2 possibilities, depending on if the numerator in (|4.4p van- 
ishes for some positive integer x or not. 

Hahn polynomials. Let (|4.4p vanish for x ~ N + 1 but be positive for x < N. 
Then one of ai and 0:2 must be equal to —N — 1, say a2 = —N — 1 and we can 
rewrite (|4.4p as 

Wjx) ^ {j: + a,)(N + l-x) 
W{x-1) x{(3i + l + N-x) ' ^ ' ^ 

where we denoted /3i = —013 — — 1. Thus, up to a constant factor, the weight 
becomes 

^(-^=[ , )[ N-x 
The corresponding orthogonal polynomials on — {0,1,..., iV} are the Hahn 
polynomials Qn{x] ai, /3i, N) given by 

Qn{x;ai,l3i,N)^3F2i ' ^^^^^^ (4-6) 
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In order to have a positive weight at aU x G the right-hand side of (14. 5p must be 
positive for a; = 1, 2, . . . , TV. Putting x — 1 and x — N we see that (ai + l)(/3i+iV) > 
and (ai + N)(l3i + 1) > 0. Starting with /3i > -1 or /3i < -1 shows that there 
are two sohitions, 

(i) ai > -1, /3i > -1, 

(ii) ai < -iV, Pi < ~N. 

Conversely, it is clear that if ai and /3i satisfy (i) or (ii) the weight function W{x) 
will have a constant sign on . Hahn polynomials satisfy the following orthogonal 
relation 



^ (ai + l),(/3i + l)jv- 

x=0 



x\{N-xy. -^nix;ai,puN)Qmix;ai,l3i.,N) (4.7) 

_ (-l)"n!(/3i + l)n{n + ai + A + l)jv+i 
iV!(2n + ai+/3i + l)(-iV)„(ai + l)„ "'™" 
The difference operator D takes the form 

D = x{Pi +N + 1- x)AV + {{ai + l)N - x{ai + /3i + 2))A. (4.8) 

Hahn-type polynomials on Nq. Assume now that the right hand side of (|4.4p 
does not vanish for any a; G N. Then it must be positive for all x G N. In order 
to have a positive function W(x) and to have convergent series X^aieNo ^^^^{^) for 
some n G No, the parameter 0:3 must be positive and ai, a2 must satisfy one of the 
following 

(i) ai and 012 are both positive, or there exists a negative integer k such that 
ai, ^2 e (k, k + 1); 

(ii) a2 = ai G C \ M. 

For such ai, a2, as we can write the weight W{x) as 

m x!(a3 + l). • ^^-"^ 

The series X]i;>o^"^(^) converges absolutely for n < aa — ai — Q!2 — 1, i.e. the 
corresponding polynomials will be orthogonal on No up to a given degree. It seems 
that these polynomials have not appeared in the literature before. Therefore, we 
derive below their basic properties: explicit formula in terms of hypergeometric 
functions, the three term recurrence formula, etc. 



Proposition 4.1. The polynomials 

1 - 



u t \ I? I n,n-a3 + ai+a2 + l,-x , /a Tn\ 

i?„(x)=3F2 + l ;1 ' (4-10) 



satisfy 

D{Rn{x)) = n(Q!3 - ai - Q!2 - 1 - n)Rn{x) (4-11) 

where 

D = -x(x + a3)VA + ((^3 -ai-a2- 2)x - (1 + ai)(l + a2))A, (4.12) 
and the orthogonal relation 

f: i^„(.)E^(.) ("^ +^^- =0, 

for n ^ m such that n + m < as — ai — a2 — 1 ■ 
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Proof. We look for a solution of ()4.1ip of the form 



i?„(x) = ^afcTOfe(a;), where nikix) = . (4-13) 

k=0 

Notice that 

Am/j = — m/j_i, xVnik — kruk, and xnik-i = {k — l)mk-i — kruk- (4-14) 
Using these formulas we can easily calculate D(Rn): 



D{Rn{x)) + a^)x\I {akmk-i{x)) - ((as - ai - a2 - 2)x 

k=l 

- (1 + ai){l + a2))akmk-i{x)] 

n 

= '^[{k - a3 + ai + a2 + l)x + a3(fc - 1) + (1 + ai)(l + a2)] akmk-i{x) 

k=l 
n 

= ^ [k{a3 — ai - Q!2 - 1 — k)ak + {ai + k + l){a2 + k + l)afe+i] , 

k=0 

where in the last formula we assumed that a„+i — 0. Thus the difference equation 
(|4.1ip leads to the recursive relation 

_ (n — fc)(n + - as + ai + a2 + 1) 
" (ai + fc+l)(a2 + fc + l) 

whose solution with aq = 1 is 

(-n)fc(n - as + ai + a2 + l)fc /.-,r\ 

^ (ai + l).(«2 + l). ■ ^'-''^ 

This shows that the polynomials Rn{x) defined by (|4.10p satisfy (|4.1ip . The or- 
thogonality follows from Remark 13.41 □ 



From the explicit formula we can derive other properties of these orthogonal 
polynomials. 

Proposition 4.2. The polynomials i?„ satisfy the three-term relation 

xRn{x) = a„i?„+i(a;) - (a„ + c„)i?„(a;) + c„i?„_i(a;), (4.16) 
where the coefficients are given by 

{n + ai + l){n + a2 + !)("- - as + ai + a2 + 1) /, x 

On = 7 TT r (4.i/aj 

(2n — as + ai + a2 + l){2n - as + ai + a2 + 2) 



_ n(n - as + ai)(n — as + a2) 

(2n - as + ai + a2)(2n - as + ai + a2 + 1) 

Proof. Using again (|4.14p we see that 



(4.17b) 



xRn = ^k{ak - ak-i)mk{x) - {n + l)a„m„+i(a;), 

k=l 

where Ofc are given by (|4.15p . It is now a simple matter to compare it with 
a„i?„+i(a;) — (a„ + c„)i?„(a;) + c„i?„_i(a;) and verify that the coefficients of mk(x) 
agree for fc = 0, □ 
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Proposition 4.3. The norm of the polynomial Rn can he computed from the for- 
mula 

nl 

{Rn,Rn) (Q,g _ Q,^ _ - 2n - 1)(1 + ai)„(l + a2)„ 
^ r(Q;3 - Qi - a2 - n)T{a3 + 1) 
r(a3 — ai — rt)r(Q!3 — a2 ~ n) 

Proof. Using the three-term relation, a standard computation shows that 

{Rn,Rn) — !!-ILJ: — (^Jl^ ^ J{g'^ ^ 

ari-ia„-2 • • ■ ao 

The quantity (Rq^Rq) = (1, 1) is a 2F1 evaluated at 1, which can be computed in 
closed form using Gauss' formula 

-j^^ _ r(Q;3 + l)r(a3 - ai - a;2 - 1) 
r(a3 - ai)r(Q!3 - Q!2) 
The statement follows from the above relations, (I4.17P and the formula 
r(a;) = (-l)"(-a; + l)„r(x-n). 

□ 

Remark 4.4. Notice that although (i?„,i?,„) converges only for finitely many 
(n, m), formula (j4.10p defines polynomials for every n g Nq. The operator D given 
in (|4.12p is admissible according to Definition 13.31 Moreover, equations (|4.1ip and 
(|4.16p hold for every n £ No- In particular, this implies that the polynomials 
Rn{x) provide a solution to the discrete-discrete version of the bispectral problem 
discussed in [3j. 



4.2. Linear eigenvalue: a = and 5=— 1. Ifa = and 6 = — 1 equations ()4.2|) 
give A{x) = Ix, B{x) = -x + s, and j3{x) = A{x) + B{x) = (/ - l)a; + s. Thus 

W{x) ^ {l-l){x-l) + s 
W{x-l) Ix ■ ^ ' 

We have now 2 possibilities: / = 1 or / ^ 1. 

Charlier polynomials. If Z = 1 equation (|4.19p essentially reduces to 

W{x) _ s_ 
W{x - 1) ~ x' 

i.e. s must be a positive number and (up to a constant factor) W(x) — /x\. The 
corresponding polynomials are the Charlier polynomials 

-n, ~x ^ 1 



Cn{x;s) = 2F0 
orthogonal on Nq. The difference operator takes the form 



D = xAV + (s - x)A = sE-{x + s)I + xE^K (4.20) 
Assume now that Z / 1. Then, we can rewrite (|4.19p as 

=c^, (4.21) 



W{x-1) X 

where we have put c = {I — l)/l ^ and d ^ {1 + s — l)/{l — 1 
We now have two possibilities depending on the sign of c. 
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Krawtchouk polynomials. If c < then (|4.21|) can be rewritten as 



W{x - 1) 

from which it follows that (— d) must be a positive integer (otherwise, W{x) will 
become negative at some point). If we put d = —N — 1 and p = c/(c — 1) G (0, 1) 
we obtain the Krawtchouk polynomials 



orthogonal on — {0,1,..., A^} with respect to the weight W{x) — p^{\ — 
p)^^^(^). The difference operator takes the form 

D^x{l-p)AV + {Np~x)A. (4.23) 

Meixner polynomials. If c > then W{1) > implies that P := d + 1 > 0. 
In order to have convergent series, we need c < 1. This leads to the Meixner 
polynomials 

M„(a:; p, c) = 2F1 (^""^""^ ' ^ " ^) ' ^^■'^^^ 

depending on two parameters < c < 1 and /3 > 0, orthogonal on Nq with respect 
to the weight W{x) = c^{(3)xlx\. The difference operator is 

£> = ^— AV + f-x - A. (4.25) 

1 — c \ 

The Meixner polynomials satisfy the orthogonal relation 



^ i^c-M„(x,/3,c)M„(a:,/?,c) = (4.26) 



c "nl 



5. MULTIVARIABLE CASE WITH QUADRATIC EIGENVALUE: a = -1 

We first solve the compatibility conditions (I2.12p in arbitrary dimension d > 1, 
and then we write explicitly all possible orthogonal polynomials with the corre- 
sponding weights. For a = — 1 formulas (|3.16p - (|3.18p give: 



"ij = Xj{xi ~ /ij) (5.1a) 

d d 

Pi^ ~Xi'^Xk + '^li,kXk+bxi + St (5.1b) 



' d d \ 

- ^ Xfc + ^ lk,^ 
V, k=l k=l J 



li = Xi \ - ^_^Xk + ■ (5.1c) 



The next proposition gives the general solution to the compatibility conditions 
Proposition 5.1. For a = —1 there are 2 solutions to the compatibility conditions. 
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Solution 1: 



d 



Pi = -{xt -h) iy^Xk~b~r 
\k=i / 

(d d \ 

k=l k=l / 



li = -X 

where b, r, and {li}f^i are free parameters. 
Solution 2: 



(5.2a) 
(5.2b) 

(5.2c) 



— Xj(^Xi li^ 

/ d d 

f3i ^ ~{xi - k) y^xfc - y^^fc + 1 

\k=l k=l 

Cd d \ 

- ^Ik - n j , 
fe=l fe=l / 



(5.3a) 
(5.3b) 

(5.3c) 



where {li, ri}f^i are free parameters and b — ^fe ~ 1- 



Proof. Let us take i j and let k ^ i,j. Comparing the coefficients of x^XjXk and 



and 



on both sides of p.l2ap we get: 



{lj.k + lj., + 2) = -2{lj., + l), 



\k=l 



ik^ + 1) 



J2 h.j ]-ib + k,, + 2) 



\k=l 



(5.4) 



(5.5) 



equation l|5.4p simply means that for i ^ j, Ij^i is independent of i, i.e. we can put 
Ij :— Ij^i for j ^ i. Using this in (|5.5p we obtain that 

= h{h -h-3 (5-6) 
Our formulas simplify as follows 

a.ij = Xj(xi - k) (5.7a) 

d 

f3t = -Xi ^ Xfe + ^ Xfc + (Zi,i + b)xt - li{li^i + b- k) 



k=l 



k^i 



/ . 

Xi \ -^Xk + ^lk + k,i 
\ k=l k^i 



(5.7b) 
(5.7c) 
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With the above formulas it is very easy to solve completely (|2.12a|) . Indeed, after 
canceling the common factor XiXj{xi — h — l)(a;j — Ij — 1) equation ()2.12ap gives 

ih,^ -h- Ij.j + Ij) Ik-b-l^^Q. (5.8) 

Let us denote — li^i — li. Then the first factor in the last formula is simply — rj. 

From (15. Sp it follows that we have 2 possibilities: — rj for all i ^ j, i.e. r :— r.i 
is independent of i, or b = X!fc=i ^fc ~ 1- The first leads to Solution 1, and the second 
leads to Solution 2. It is immediate to see that the functions defined by (|5.2p or 
([Q]) satisfy ([2l2bl) . □ 



Next, we want to determine the possible weights and orthogonal polynomials. It 
is easy to see that there is no weight corresponding to Solution 2. Otherwise, from 
([53]) and l[2J|) we will have 

-"'^^^-^or^ = l,2,...,d. (5.9) 



W{x - a) Xi 

Then, for a; = it follows that 1^ < 0, i.e we can put — — Vi > Q and therefore, 
up to a constant factor, W{x) is given by 



w{x) = n 



i—l 

but then even (1,1) will diverge. 

Let us now concentrate on Solution 1, given by (|5.2p . Equations (15. 2p and p. 9 
give 

W{x) x,-k^l (ELi - 5 - r - 1 



(5.10) 

for i = 1, 2, . . . , d. For a vector y G K*^ we will denote 

fe 

i2/i = Ey- (5.11) 
1=1 

a, = -(/, + 1), (5.12) 

o- = (cri,(T2, . . . ,crrf) (5.13) 

iy(x) Xi + ai (|x|-&-r-l) 



If we put 
and 

we can write (|5.10p as 



W{x-ei) X, {\x\ + \a\+d-ry '■^'"^^^ 

We have several possible sub-cases, depending on whether the right hand side in 
equation l|5.14p can vanish or not. We divide them first in two big subclasses: 

• & + r ^ No; 

• & + r e No. 



5.1. Case 1: 6 + r ^ Nq. Then \x\ — h — r — 1 does not vanish. We have two 
possibilities: 
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5.1.1. For some i G {1, 2 . . . , d}, (— o"i) N. We show below that this impUcs that 
(Tj > —1 for aU j G {1,2..., d}. Take j ^ i, arbitrary fc G N and n G Nq. From 
(j^TTi)) we see that 

W{kei + nej) _ k + a, k + n- b- r- l ^ ^ 
W{{k -l)ei + nej) ^ k k + n+\<T\+d~r ' \ - ) 

The second ratio is positive for large n, and therefore k + ai > for every fc G N, 
i.e. Ui > —1. Equation (|5.15p for n = implies that 

k — b — r — 1 
k + \a\ + d - r 

for every k €N. But then for every j G {1, 2, . . . , d} we must have 

Wjcj) _ 1 + gj 1-b-r-l 
W{0) ^ 1 l + |(T|+d-r ^ ' 

showing that (jj > —1. 

Since aj > —1 for all j and 6 + r ^ Nq we must have V = Nq. Up to a constant, 
the weight Vl^(a;) is given by 



X,! (|cr|+d-r + l)|^| 
Changing the parameters 

-f^\a\+d~r /? -(6 + r) - 1, 
the weight takes the form 

and the free parameters are {ui}, /3 and 7. The difference operator becomes 
D = - ^ a;j[cri + 1 + + (5ij(7 - |cr| - d)]AiVj 

l<ij<d 

d (5.16) 

+ - /? - |a| - - 1) - (1 + /3)(1 + fxOlA,, 

i=l 

and A„ = n{—n + 7 — /3 — |(t| — d). 

For a vector y — (yi,y2, ■ ■ ■ ,yd) G K'^ we will denote 

= {yj,y3+i, ■■■.Vd) and y, = (2/1,2/2, • ■ • ,2/i), (5-17) 

with the convention that y'^'^^ = and Yq = 0. Let us denote by Rn{x; ai,a2, as) 
the orthogonal polynomials given by (|4.10p . 

Theorem 5.2. For v G Ng, such that 2\v\ < 7 — |cr| — — d — 1, i/ie polynomials 

d 

R^{x;a,P,j) = ]^(q!2.j + l)i,^R^^{xj;aj,a2j,a3j), 
i=i 

where a2j and a^j are given by 

a2,j=f3+W'+'\ + \Xj^i\ 

«3,, = 7 - V K+'l - (d - j) + l^.-il 
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satisfy the difference equation 
and the orthogonal relation 

_ (1 + /3)|,|r(l + 7)r(7 -f3-\a\- 2W\ - d) 

r(7 - H - IH + 1 - d)r(7 - /3) 

d u \ f , T\ (5.19) 

Proof. Notice that aa.j — a2.j is independent of x, and therefore from (|4.10p one 
can immediately see that 4>i, is indeed a polynomial of x of total degree From 
Theorem 13.71 we know that for every k S No, the equation Du = XkU has — 
C^^fc"^) ~ dim(nj!/n^_j^) linearly independent solutions. Therefore, it is enough 
to prove that (|5.18p holds. 

In (jS.lSp we will first sum with respect to Xd, then with respect to Xd-i, and so 
on. Writing 

and extracting only the terms depending on Xd we get the sum 

T --S^R (n^\R ^''d + lU {[3+l + \Xd-i\)x, 

'"i^o (7+l + l^.-i|)x. 



where Rm,{xd) — Rm{xd', crd, a2,d, ct^^d)- Using now (|4.18p and the fact that a^^d — 
(^2,d — 1 ~ P is independent of x, we see that 

J r(7 + |x,_i| + i) 

^^^^^ |X,_i|),,r(7 -ad-i^d + \Xd-i\) ' 
where B,^^ is a constant (independent of x), whose value can be extracted from 

^ i/d!r(7 - P - <7d- Vd) 

""^ ~ (1 + nd)uAl -l3-ad-2vd- l)r(7 - P ^ i^d)' 
Using the fact that 

T{x)^{x-n)nT(x~n), (5.20) 

we see that 

J _ c ^ (7+ 

(1 + /? + l^.-i|)..(7 -^d- ^d)\xU 
with = B^A^(^ + l)/r(7 — Od — Vd)- This shows, in particular, that 



(7+ 

~ "'^'^ '''^'^^^'''{l-^d-^d)\x,_,\ 



(5.21) 
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where we have used the identity 

(/3 + + 1 + \Xd-i\)., = {P + l)|.,|(/3 + 1 + yd)\x,_,\- (5.22) 

Equation (|5.2ip shows that the remaining d — 1 fold sums of (i?^ , R^) have exactly 
the same structure as that of the original d fold sums with /? and 7 replaced by 
(3 + Vd and ^ ~ Ud — Vd — respectively. In other words, it shows that we can use 
induction to complete the proof. For we obtain the following formula: 

(i + /?)i.ir(i + 7) -A uj\ 

" r(7-la|-kl + i-d) -Li (1 



r(7 - \a\ -\y\ + l-d) f \ (1 + CT,),^ (7 - /? - \ai\ - 2|W-| + j - d - I) 



^-Q T{^-f3~\a^\-2\:.^+'\-,y,+j-d) 



LJ. r(7 - /3 - |a^-+i I - 2|i.^+i I - + - d) • 

(5.23) 

Using several times (|5.20p one can rewrite the right-hand side as in (|5.19p . □ 

5.1.2. For every i G {1, 2 . . . , d} we have (-CTi) ^ k + 1 e N, i.e. k £ Nq. The 
corresponding V is the parallelepiped 

By (|5.10p the weight function in this case is given by 



d 



"-■M-n^ ,_ir.tl),., - ^''-'> 

where we again set /3 = — (6 + r) — 1. Putting x = Ci and x = {li,l2, ■ ■ ■ ,ld) hi 
([5?T0l) we get that 

0+1 \l\ + B 

^ > and ^ > 0. 



1^1 + r - 1 r 
From this we see that the parameters (3 and r must satisfy one of the following 
conditions: 

(i) /3 > -1 and r > 0; 

(ii) /3 < ~\l\ and r < -\l\ + 1. 

Recall that b + r ^ Nq, but everything will hold even if5-|-r = — /?— 1 £ No as long 
as /3 + 1^1 < (i.e. if (ii) holds), because (/9 + l)|j.| will not vanish for x G Vi'^. 

In the following we will use the notations and Yj defined in (|5.17p . For 
example, \Lj\ = li + ■ ■ ■ + Ij. If ^ Ng, then v < I means j/^ < for 1 < i < d. 
Recall that the Hahn polynomials are denoted by Qn{x; ai, /3i, N). 

Theorem 5.3. Let U E Nq, 1 < i < d. For v G Ng, Vi < li, the polynomials 

d 

cj)^{x][i,r,l) = ]^(aij + l)^.Q„^{xj;ai^j,a2,jJj), (5.25) 
1=1 

where ct\_j and 02, j are given by 

ai., + and a2,, = |L,-i | - | + | + r - 1, 

satisfy the difference equation 
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and the orthogonal relation 

^(j)^{x; f3,r,l)(j)f,{x] (3,r,l) f[ ^ ^'j'"'' ..^J^ ^''l — = B^S^^^,, (5.26) 

where the normalization constant is given by 

(-1)H(1+/3)|,| ■A ^^.!(/3 + r + 2|^^+i| + i/, + |L,_i|),,+i 
^ (^+IH)|;|-M /J H,).,(/3 + r + 2|^^| + |L,_iI) ■ ^ - ^ 

Proof. Since j + a2,j is independent of x, it is easy to see from (I4.6P that is 
indeed a polynomial of x of total degree We proceed as in the proof of Theorem 
15.21 In (|5.26p we will first sum with respect to Xd- Using the fact that 



(-1^1 - r + = (-1^1 - r + 1 + \Xd-i\h,i~\l\ - r + 
iP + 1)|.| = (/3 + + 1 + 

we can split the weight function as a product, 

_ {-ld)xAP+^+\Xd-l\)x^ Mill IV \ 

^\^) - — rnn 1 1 I iy — h — ^ i^d-i) 

a;d'(-K| - r + 1 + 

where 



It is easy to verify that 



Hence, using the fact that a\^d and a2,d are independent of a;^, the sum over Xd in 
(I5.26P becomes 

E''' ^ . . X /a:^d + ai.<i\ /^d - a:^ + Q;2,d 

QvA^d)Qp.A^d)\ , _ 

x^=o ^ ^'^ ^ \ «d a;d 

where Qm(x<j) = Qm(a;d, ai,d, a2,d, ^d)- Using (|T7)) and simplifying, we get 

r i-ir^iydljiyd + 13+ \Ld-i\ + r);,+i (aa.d + 1)., 

{~ld)vdV'^d + P + \Ld-i\ + r) [a2,d + Ijid 

In the above expression, the term (ai^^ + 1)^^ contains variables xi, . . . , Xd-i- Com- 
bining this term with W'{Xd-i) and using (|5.22p . we see that the weight function 
for a;i, . . . , Xd-i becomes 



,1 a;«! - r + l)\Xd-i\ i'^2,d + 
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By the definition of a2,d and expanding the Pochhammer symbols gives 

(-1^1 -r + l)|^^_^l ' ' 
(a2,<i + l) 



= {-iy^--'{-\l\ - r + -r + 2)... {-\Ld-i\ - r - j^d + \Xd-i\) 

= HLd-i\ - r + 1 - i^d)\Xa_,\{\Ld^i\ + r + i^d)!^-^^. 
Consequently, the weight function Wd-i becomes 

^ + "fr i-hU (/? + l + ^rf)|X,_,| 

Wd-i[^d-i) - — ; T I I j — -^-j7 j — ^ ■ 

[\Ld-i\+r + iyd)i^-ud fj^ x^\ (-\Ld-i\ - r + 1 - iyd)\x^_,\ 

Apart from a constant multiple, the weight function Wd-i has exactly the same 
structure of 14^(2:) with (3 replaced by P + Vd, f replaced hy r + Vd, I replaced by 
Ld-i, respectively, and one variable less. Hence, we can use induction to complete 
the proof. For B^, we obtain the following formula 

R-r-nMn+m n ^,!(/3 + r + 2|.^-+i|+., + |L,-_i|),^.+i 

which combined with 

d 

Y[i\Lj-i\ + r+ = (r+ 

gives (ESZl)- □ 

5.2. Case 2: h + r = N ^ Nq. Again we have two possibilities: 

5.2.1. For every i e {1, 2 . . . . d}, li + l ^ {1,2,..., TV}. In this case the numerator 
in (|5.14p is zero when |a::| — N + 1 and we have V — V^. The corresponding 
orthogonal polynomials are the Hahn polynomials studied in _4_ . Setting Ui = — 1 
for 1 < z < d, cr = ((Ti, .. . ,0-d+i), 

r = iV+|cr|+rf+l, and 6 = -(|cr| + + 1), (5.28) 

equation (|5.14p gives 

W{x) _ x,+a, N + 1- \x\ 
W{x-e^) ~ X, N + l + Gd+i-\x\' 

The above ratio must be positive for all x G V^. In particular, for x — e-i and 
X = Nci we see that 

"'^^ >0 and ^1±J^>0 fori = l,2...,d 



N + ad+i 1 + (Jd+i 

From this it follows easily that the parameters {o'i}'!^^^ satisfy one of the following 
conditions: 

(i) CTi > -1 for i = 1, 2, d + 1; 

(ii) cr, < -N for i 1, 2, . . . , d + 1. 

The weight function takes the form 

Xi + (Ti\ /N - \x\+ ad+] 
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with {(Ti} as free parameters. If (T; > — 1 for i = 1, 2, . . . , + 1 or if (Ti < —N 
for i = 1, 2, . . . , d + 1 but iV is even then W{x) > on V^. If a, < -N for 
i = 1,2, . . . ,d + 1 and if N is odd then W{x) < on V^, so one needs to change 
the sign in formula (|5.29p . in order to get a positive function. 
The difference operator takes the form 

d 

D = "^x^{N -Xi + \a\~a, + d)A,Vi- ^ Xj{x, + ai + l)AiVj 

d 

+ Yy^^ - ^^(C^^ + 1) - X^(\(J\ - (T; + d)]Ai 
i=l 

and the eigenvalues are A„ = —n{n + \a\ + d). 
Theorem 5.4. For v eW^ and \v\ < N, the polynomials 

Q,{x;a,N) = t^T\ ^1l±^^^n + \X,^, \ + \v^+^\),^ (5.30) 



(-^)m /J 

xQ.,(x,;a,,a„iV-|X,_i|-|i.^+i|), 
where Oj — \a^^'^ \ + 2|i^-'+"'^| + d — j, satisfy the difference equation 

and the orthogonal relation 

^ QUx;a,N)Q,{x;a,N)f[('''^''')(^ -^''l+J^^^^ (5.31) 



|2;|<JV i=l ^ 

where Ay is given by 

d 



(5.32) 



These formulas are essentially contained in [3] . They can be deduced from (|4.7p 
as in the proof of Theorem 15.21 Explicit biorthogonal (not mutually orthogonal) 
Hahn polynomials were also found in [8]. 

5.2.2. There is a nonempty set S C {1, 2, . . . , d} and k + I G {1, 2, . . . , A^} for 
i & S. In this case the set V is 

V^^S ^V^n{x: X., < h for i G S}. (5.33) 

If 5 = {1, 2, . . . , c?} we can also assume that h + ■ ■ ■ + Id > N, otherwise it becomes 
the parallelepiped case discussed in Theorem 15.31 

The weight W{x) is again given by (lEMl), but is replaced by V^^. The 
corresponding polynomials are the same as in Theorem 15.41 with the restriction 

G V^^s- 

Theorem 5.5. For v G 5 the polynomials Q,^{x;a,N) defined by (I5.30p satisfy 
the difference equation 
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Xi + at 



N - \x\ 
N- 



(Jd+i 



and the orthogonal relation 

d 

Y,Q,{x-a,N)Q^{x-<j,N)\{ 

x£V i=l 

where Ai, is given by (|5.32p . 

Remark 5.6. Notice that (|5.25p in the paraUelepiped case V^'^ gives a polynomial 
of total degree We can use the generators Gf in Example 13. II to write the same 
polynomial using only the monomials x^ with /i G Vf". Similarly, in the case of 
Theorem 15.51 we can use the generators {{—Xi)i-j^i : i Cz S} to express Qi, in terms 
of the monomials a:^ with /i G g. 

Remark 5.7. In Theorems [Oll^ we have = |Afe(y)| < dim(n^/n^_J for 
some fc's and therefore the equation Du(x) = \u{x) will hold a priori only if we 
consider u{x) as an element of However, one can show that the corresponding 

polynomials can be obtained as a limit from the polynomials in Theorem 15. 21 This 
fact can be used to show that the equation Du(x) = \u{x) actually holds in the 
space R[a;i , a;2 , . . . , Xd\ ■ 

5.3. Summary. If a = —1 we have essentially one difference equation with coeffi- 
cients given in (|5.2p . By specifying the free parameters, however, we end up with 
four different types of solutions given in Theorems 15.21 - [531 respectively. 

As an example, let us consider the case d = 2. The weight functions and the 
corresponding sets on which they live are listed below: 



(i) W{x) = 

(h) W{x) = 

(in) W{x) = 

(iv) W{x) = 

(v) Wix) = 



Xi\x2\ (-|cr| 



V 



Xi\x2\ (-/-r)|j,|' 
a-i + ^)xi{o-2 + 1)0:2 (^3 + ^)n-\'x\ 



+ 1) 



Xi\x2^- 



i-hh, (f^3 + 1) 



iN-\x\y. 



V = v, 



2 . 

N' 



N~ 



xi\ X2\ {N-\x\y. 

-h)xi (0-2 + 1)2:2 ('^3 + l)Ar-|a:| 



V = V^ s, where S = {!}. 



xi\ X2\ {N-\x\)\ 
In the last case, one can also exchange xi and X2 to get another case. It should be 
mentioned that the equation (|5.2p was considered in [T^, but only the case (iii) 
was identified there. 



6. MULTIVARIABLE CASE WITH LINEAR EIGENVALUE: a = AND 6 = -1 

When a = and 6 = — 1, equations (|3.15p - (|3.18p become 

a^j = -h,jXj for i / j (6.1a) 
Bi^-Xi + Si, (6.1b) 

A = ^ li,kXk + {hs - l)xi + Si (6.1c) 
d 

li^Xi^lk,i- (6. Id) 
fc=i 
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First, notice that we cannot have 7^ ~ 0. Indeed, if we assume that 7^ = then 
in order to have a self-adjoint operator we need (3i = 0, which imphes that h^k = 
for k ^ i. But then ai^k = for k ^ i and therefore ak.i = for k ^ i, which 
simply means that the operator D is independent of Xi and . Thus we have an 
operator acting in a (d — 1) dimensional space, trivially extended, by adding Xi as 
explained in Remark 13.61 

Below we assume that 7^ 7^ for all i. Notice that 7^ depends only on x^, so in 
(j2.12bp we can cancel 7^7^ and we get 

I3i(x - ei)(3j{x - Ci - Cj) = (3j{x - ej)(3i{x - - Cj). (6.2) 

But Pi{x — Ci) = (3i{x) + 1 — and (3i{x — Ci — Cj) — f3i{x) + 1 — li^i — Uj and 
plugging these in (|6.2p we get 

,(A(x) + 1 - k^i) = k,j{Pj{x) + 1 - Ij^j), (6.3) 

for i ^ j- Comparing the coefficients of Xk in the last formula we see that 

l].ik,k = k,]l].k for k ^ij (6.4) 

These two equations have essentially two different types of solutions, which we 
discuss in two subsections. 

6.1. Case 1: 1,^^ ^ ior all 1 < i j < d. 

Proposition 6.1. Assume that ^ for all 1 < i ^ j < d. Then the most 
general solution of the compatibility conditions ()2.12|) is given by 

aij — —liXj for i ^ j (6.6a) 

Bi — ~Xi + Us, (6.6b) 

d 

|3^=l^\yXk+s\ (6.6c) 



(6.6d) 



Proof. From equation (|6.5|1 it follows that = — 1. Denote Ij :— Ijj — 1. 
Then we have = Ij for all i ^ j- Using now (|6.3p we see that IjSi = Usj. Thus 
Si = hs, which leads to formulas (|6.6p . 

Conversely, it is straightforward to see that if we define aij , Pi and 7^ as in (|6.6ap , 
(|6.6cp and (j6.6dp . then the compatibility conditions in Corollarv 12.31 are satisfied, 
i.e. the above formulas give the most general solution in the case lij 7^ 0. □ 

Below we determine the weight functions and the corresponding orthogonal poly- 
nomials. For every i we have 

W{x) k \x\+s-l 




where s and {li}f^i are free parameters. 



W(x^eC) 1^1 + 1 

and therefore 

W{e,) hs 



>0, 



1^(0) \l\ + l 

which shows that li must have the same signs. There are two possible cases. 
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6.1.1. k — —pi < for all i. In this case, 

W{x) Pi 1 — s — \x\ 



W{x ~ Ci) 1 - \p\ X., 

If the denominator does not vanish, then for large the second ratio will be 
negative and therefore we must have 1 — |pl < 0, i.e. \p\ > 1. But then if we denote 
Ci = Pi/{\p\ — 1) we will have |c| > 1 and up to a constant factor the weight is 
W{x) — (s)|:r| HiLi c^V^^jii which leads to divergent series. 

Hence the only possibility here is (— s) — N ^ No. This forces \p\ < 1 and 
V = V^. The difference operator is then 

d 

D= J2 -ft)a;jA,V, +^(p,iV-a^,)A«, (6.7) 

l<ij<d i=l 

and the eigenvalues are A,i ~ —n. The orthogonal polynomials are the Krawtchouk 
polynomials on . Recall that Krawtchouk polynomial in one variable is denoted 
hyKn{x;p, N). 

Theorem 6.2. Let < p, < 1, 1 < i < d, and \p\ < I. For v S N[J, \v\ < N, the 
polynomials 

K^{x;p,N) ^j^^ {[ J. ^ —{-N+\X,^,\ + (6.8) 

satisfy the difference equation 
and the orthogonal relation 

-Pf_ (-i)i'i A "Afl 



(6.9) 



where xa+i = N - \x\, pd+i = 1 - |p| and Vd+i = 0. 

In fact, these orthogonal polynomials can be considered as a limit of the Hahn 
polynomials (|5.30p (see [SKTU]). Indeed, using the well known relation 

lim Q„(a;; pt, (1 - p)t, N) = K„{x] p, N) 

t — ^CJO 

in one variable, it is not hard to see that 

lim Q^{x;pit, . . . ,pdt, (1 -pi Pd)t,N) = K^{x;p,N). 

t — >oo 

The orthogonality of K^{x;p, N) follows from (|5.3ip under the limit. We refer the 
reader to [5j for other properties and applications of Krawtchouk polynomials. 

6.1.2. k > for all i. In this case s > 0. Denote Ci = li/{l + \l\). Then < 1 and 
|c| < 1, the weight function is 

mx) = ^c^ = (,s),.,nj:. 

2—1 
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and y = Nq. The difference operator takes the form 

^= E (^»^. + t^)^.^'V. + E(--' + t^-)'^'- (6-10) 

The orthogonal polynomials are the Meixner polynomials on Nq but they are dif- 
ferent from product Meixner polynomials. Recall that Meixner polynomial in one 
variable is denoted by M„(a:; (3, c). We write \Cj \ — Cj + c^+i + ■ ■ ■ + Cd and define 
Cd+i = 0. 

Theorem 6.3. Let < Ci < 1, 1 < i < d, and |c| < 1. For v G Ng, the polynomials 

where Sj = s + {v-'^^l + \Xj_i\, satisfy the difference equation 
and the orthogonal relation 



(6.12) 

Proof. These relations can be derived inductively using the orthogonality (j4.26p of 
the Meixner polynomials of one variable. It is clear from ()4.24p that defined by 
(|6.1ip are indeed polynomials of x. We can write the weight function as a product 

W{x) = ^l!li±fel^T4^'(Xrf_i) 

where 

w^'(x,_i)=(.s)ix._.inf^- 

Using (|4.26|) the sum over Xd in (|6.12p becomes 

,cT{s + \Xd^i\)., 



= + \Xd-i\),,{s + \Xd-l\)^,, M,A^d)M^^{ 



n ^d- 

Xd=0 



''''''' {I - CdY (l-cd)l^-il - 



where M™(xrf) = M„(a;d,s + |Xd_i|,Cd). 

Combining this with W' we obtain the new weight function 



1=1 



Apart from a constant multiple, the weight function Wd-i has exactly the same 
structure of W{x) with s replaced by s + Vd, and Ci replaced by Ci/(1 — Cd) for 
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i — 1,2, . . . ,d — 1. The proof now follows by induction. The constant {M„, M„) is 
given by 



which leads to (|6TT2| . □ 

Remark 6.4. It is well known that the Meixner polynomials are the limit of the 
Hahn polynomials, 

1 - c 

lim Qn{x;h-l,N , iV) = M„(x; &, c), 

N^oo C 

There is an analogous relation between the Meixner polynomials (|6.1ip of several 
variables and the orthogonal polynomials on the parallelepiped in (|5.25p . at least 
when parameters c have rational values. Indeed, let (j)u{x] f3,r,l) denote the poly- 
nomials defined in (15.25^ . which are orthogonal with respect to the weight function 
in ()5.24|) . In these polynomials we set /? = s — 1 and 

Using induction, it is easy to see that Ij/cj — Ij-i/cj^i, which then implies that 



\Lj-i\ +r= [ ^ Ij, l<j<d. 



Ij ^i\l\+r)c,, or l^^^-lj-^r, I < J < d. 

If Cj are rational numbers, we can choose r = {1 — \c\)N, N G No, for certain N, 
so that Ij are integers. Upon taking r — > oo, which shows Ij oo, it follows from 
UjAM^i that 

lim (f>^{x; s — 1, r, /) = M^{x; s, c). (6.13) 

I — >oc 

6.2. Case 2: = for some I < i j < d. Recall that aij = for i ^ j \i and 

only if aj^i = 0. This essentially means that and Ij^i are simultaneously zero or 
nonzero for i ^ j. 

Lemma 6.5. If j k and Ij^k — then for every i ^ j,k we have li j = Ij i = or 

Proof. Follows immediately from (|6.4p . □ 

Theorem 6.6. Assume that = for some 1 < i j < d. Define I — {i} U {m : 
lm,i 7^ 0} and J = {1, 2, . . . , 0?} \ / D {j}. Then D = Dj + Dj, where Dj is an 
admissible operator in the variables {xm : m G /} and Dj is an admissible operator 
in {xk : k ^ J}. 

Proof. First we show that for m d I and k ^ J we have lm,k — Ik.m = 0. Indeed 
if m = i, then Ik^i = by the definition of J. If to ^ i then Im^i ^ 0, but Ik.i — 0. 
The previous lemma shows that lk,m = lm,k — 0. 
Thus if TO G / and fc G J we see at once that: 

• Ctm,k = ak,m = 

• 13m contains only the variables from / 

• (3k contains only the variables from J. 

The decomposition of D follows immediately from the above observations. □ 
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Remark 6.7. Theorem 16.61 savs that if kj — for some i ^ j, then the operator 
D sphts as a sum of 2 operators of independent variables. Conversely let Dj be an 
admissible operator in the variables x' = {xm : m G /} and Dj be an admissible 
operator in x" = {xk : k £ J} with / n J = and / U J = {1,2, ... ,d}. Let us 
denote hy pj^{x') and p'l{x") the polynomials satisfying 

Diipi) = ~\^'\pi and Djipl) = ~\iy\pi. 

Then if put D = Di + Dj and Pf^^^ix) = p'^{x')pi{x"), where x = {x',x") £ R'*, 
we have 

Dp^,,v{x) = + \v\)p^^^{x). 

Thus, in this case, the eigenfunctions of the difference operator are product 
of orthogonal polynomials of fewer variables, which satisfy difference equations of 
lower dimension with linear eigenvalues. For d = 1 these are the polynomials of 
Charlier, Krawtchouk and Meixner given in Section l4?2l In higher dimensions, there 
are also the Krawtchouk polynomials of several variables given in Theorem 16.21 and 
the Meixner polynomials of several variables given in Theorem 16.31 

Clearly there are many product polynomials of this type and the number in- 
creases drastically as the dimension grows. As an example, we list all cases for 
d — i below. To list the different types we use the abbreviation of C, K, M for 
Charlier, Meixner, and Krawtchouk polynomial of one variable, respectively, and 
use K2 and M2 to denote the Krawtchouk polynomials of two variables on and 
Meixner polynomials of two variables on Ng. A product polynomial is denoted by 
its components. As an example, CCM stands for a product of the type Charlier- 
Charlier-Mcixner. 

Example 6.8 {d ~ 2>). There are sixteen product types, which we list according 
to their domains of orthogonality: 

(1) nl: CCC, MMM, CCM, CMM, CM2, MM2. 

(2) X V^: CCK, MMK, CMK, KM2. 

(3) No X X 1/^, CKK, MKK. 

(4) xy^ xF^: KKK. 

(5) X No: K2C, K2M. 

(6) vk'^yk- K2K. 

Acknowledgments. It is our pleasure to thank a referee for several thoughtful 
suggestions and for bringing the reference [5] to our attention. 
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